New Constraints on the Early Expansion History 
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Cosmic microwave background measurements have pushed to higher resolution, lower noise, and 
more sky coverage. These data enable a unique test of the early universe's expansion rate and 
constituents such as effective number of relativistic degrees of freedom and dark energy. Using 
the most recent data from Planck and WMAP9, we constrain the expansion history in a model 
independent manner from today back to redshift z = 10 5 . The Hubble parameter is mapped to a 
few percent precision, limiting early dark energy and extra relativistic degrees of freedom within 
a model independent approach to 2-16% and 0.71 equivalent neutrino species respectively (95% 
CL). Within dark radiation, barotropic aether, and Doran- Robbers models, the early dark energy 
constraints are 3.3%, 1.9%, 1.2% respectively. 



I. INTRODUCTION 



II. COSMIC HISTORY MAPPING 
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Except for the last e-fold of cosmic expansion, our 
knowledge of the state of the universe arises directly only 
through measurements of the cosmic microwave back- 
ground (CMB) radiation or indirectly (as in models of 
its influence on growth of large scale structure) . Recent 
CMB data [1, 2] provides a clear window on an addi- 
tional 10 e-folds of history (back to redshift z = 10 5 ), a 
vast improvement in mapping the universe. 

The expansion rate, or Hubble parameter, is a fun- 
damental characterization of our universe, and includes 
information on its matter and energy components, their 
evolution, and the overall curvature of spacetime. More- 
over, the CMB encodes linear perturbations in the pho- 
tons and the gravitational potentials they experience, 
providing sensitivity to the microphysics of components, 
e.g. their sound speed. 

These observations lead to constraints on quantities 
such as early dark energy and extra neutrino species or 
other relativistic degrees of freedom. However most anal- 
yses assume a specific model for these deviations, en- 
abling stringent but model dependent constraints. Our 
approach is to map the cosmic state and history in as 
model independent fashion as practical, guided by the 
data. We utilize the results of the principal component 
analysis of [3] to define localized bins of Hubble param- 
eter in log scale factor that are most sensitive to the 
data, and then carry out a Markov Chain Monte Carlo 
(MCMC) analysis to constrain them. For data we use 
the most recent CMB results from the Planck satellite 
[1] and WMAP satellite [2]. 

In Sec. II we describe the methodology for mapping the 
Hubble parameter, then in Sec. Ill present the results for 
the MCMC fits of the standard cosmological plus expan- 
sion history parameters. We discuss the implications for 
dark energy, relativistic degrees of freedom, and space- 
time curvature in Sec. IV. 



For robust, model independent results we adopt a com- 
bination of principal component analysis (PCA) and bin- 
ning. This avoids assuming a specific functional form for 
the Hubble parameter or dark components and allows 
the data itself to inform where the greatest sensitivity 
lies. Such PCA on the Hubble parameter for projected 
mock CMB data was used in [3] to predict the strength 
of constraints at various epochs of cosmic history. 

This identification of the key epochs where physical 
conditions most affect the observations enables informed 
choice of bins in log scale factor to use in a MCMC fit. 
Bins have several advantages over the raw PCA: 1) they 
are localized and can be clearly interpreted physically - 
the Hubble parameter during a specific epoch, 2) they 
avoid negative oscillations that can cause unphysical re- 
sults (while the sum of all PCs will give a positive, phys- 
ical Hubble parameter squared, this is not guaranteed 
for a subset), and 3) they are well defined, not changing 
when new data is added. 

The Hubble parameter, or logarithmic derivative of the 
scale factor, H = dhia/dt, is then written as 

RttC 

H 2 (a) = — [p m (a) + p r {a) + p A ] [1 + 6(a)] , (1) 

where S accounts for any variation from ACDM cosmo- 
logical constant plus cold dark matter expansion history, 
and p m is the matter density, p r the radiation density, 
and p\ the cosmological constant density. The bins in 
the deviation 6(a) are slightly smoothed for numerical 
tractability, with 
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(2) 
Within bin i, 6 = Si and far from any bin 8 — 0. A 
smoothing length r = 0.08 was adopted after numerical 



convergence tests. (A similar binned approach was used 
in [4, 5] to bound early cosmic acceleration.) 

We modify CAMB [6] to solve the Boltzmann equa- 
tions for the photon perturbations in this cosmology. The 
dark energy density contributed by the deviations S and 
the cosmological constant term (which becomes negligi- 
ble at high redshift) has an effective equation of state 



w 



Q8 



1 + S(l + Q) 



(1 + w bg ) 



dS 



31 + S{l + Q)dlna ' 
(3) 
where Q = (p m + Pr)/pA and Wb g is the background equa- 
tion of state of the combined matter and radiation (e.g. 
1/3 during radiation domination, transitioning to dur- 
ing matter domination). Thus w and w' — dw/dlna, 
entering into the Boltzmann equations, are defined fully 
by Eq. (2) for S. We choose the associated sound speed 
to be the speed of light, as in quintessence dark energy, 
but explore variations of this later. 

Guided by the PC A of [3], we choose bins <5i_ 5 in 
the logarithmic scale factor log a = [—5,-4], [—4,-3.6], 
[-3.6,-3.2], [-3.2,-2.8], [-2.8,0] so the finest binning 
is near CMB recombination at a w 10~ 3 . The cosmo- 
logical parameters we fit for are the six standard ones: 
physical baryon density f2&/i 2 , physical cold dark matter 
density fl c h 2 , acoustic peak angular scale 0, primordial 
scalar perturbation index n s , primordial scalar amplitude 
ln(10 10 J 4 s ), and optical depth r, plus the five new devia- 
tion parameters <5i_5. Additional astrophysical parame- 
ters enter from the data, as discussed in the next section. 



III. CONSTRAINTS 

To constrain the cosmology with the data we use 
MCMC analysis, modifying CosmoMC [7]. The like- 
lihood involves the temperature power spectrum from 
the two satellite experiments, and the E-mode polar- 
ization and TE cross spectrum from WMAP (the first 
Planck likelihood release does not include polarization, 
or the high multipole likelihoods from Atacama Cosmol- 
ogy Telescope [8] or South Pole Telescope [9]; in the fu- 
ture such data should become available). Astrophysical 
nuisance parameters characterizing foregrounds (see [1]) 
are marginalized over. 

Figure 1 shows the constraints on the standard cosmo- 
logical parameters, in the ACDM case (fixing Si = 0) and 
when allowing variations in the expansion history (fitting 
for the Si). Here the Hubble constant H replaces the 9 
parameter and we omit showing r. Including the fitting 
for expansion history deviations induces roughly a factor 
of two larger marginalized estimation uncertainties for 
most of the standard cosmology parameters, and signif- 
icantly shifts the cold dark matter density value. This 
is due to the deviations in the Hubble parameter having 
similar effects on the expansion near recombination to 
those in matter, so S takes the place of some of p m . We 
discuss this degeneracy further later. The best fit for the 



ACDM case remains within the 68% confidence contour 
when allowing expansion deviations. 
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FIG. 1. Joint 68% confidence contours on the standard cos- 
mological parameters are shown when allowing for expansion 
history deviations from ACDM (black), and fixing to ACDM 
(smaller contours or dashed curves). Plots on the diagonal 
give the ID marginalized probability distributions. 

Figure 2 shows the constraints on the expansion history 
deviations. Note that to ensure positive energy density 
(and Hubble parameter squared) we restrict S > 0, i.e. 
equal or more early energy density than in the ACDM 
case (which has il\ ~ 10~ 9 at a = 10~ 3 ). Note that these 
binned deviations do not have appreciable covariances 
with each other, with the correlation coefficients under 
0.26 except for S2-S3 at 0.49. This is a useful feature 
adding near independence to localization, making their 
interpretation transparent, and is a result of the careful 
choice of bins based on the PCA of [3]. 
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FIG. 2. Joint 68% and 95% confidence contours on the expan- 
sion deviation parameters are shown. Plots on the diagonal 
give the ID marginalized probability distributions. 

Table I gives the 95% confidence upper limits on 
each expansion deviation parameter, showing that recent 
CMB data provides 2-16% constraints on the expansion 
history back to z = 10 5 . The earliest bin, of S\, is reason- 
ably constrained despite being well before recombination, 



and should improve further when adding high resolution 
(high multipole I) measurements. The second bin has 
equivalent constraints when taking into account its nar- 
rowness. Around recombination, however, S3 and 64 have 
looser bounds because all the standard cosmological pa- 
rameters also enter strongly at this epoch, and so the 
increased covariance dilutes their estimation. They have 
the two highest correlation coefficients, of 0.89 between 
S3 and Q c h 2 and —0.76 between S4 and ilbh 2 . Finally, 
the late, broad bin of 65 has strong constraints. These 
behaviors are all consistent with the pre-Planck, Fisher 
matrix predictions of [3] (see their Fig. 4). Adding late 
time data or priors (which we avoid; see concluding sec- 
tion) can shrink some uncertainties by up to 60%. 

The expansion history docs not completely define the 
system of Boltzmann equations: the effective dark com- 
ponent can have internal degrees of freedom such as 
sound speed c s that determine the behavior of its per- 
turbations and hence the gravitational clustering of the 
photons [10]. Therefore we also show in Table I the con- 
straints when this sound speed is equal to that of a rel- 
ativistic species (cf = 1/3), or is much smaller than the 
speed of light, cold dark energy with c s = 0. The c s = 
case has looser bounds, due to the additional influence on 
the photon clustering with the strengthened gravitational 
potentials, and covariance with matter parameters dur- 
ing matter domination. For the c 2 = 1/3 case, where the 
extra expansion rate corresponds to extra relativistic de- 
grees of freedom, the constraints are weaker during radi- 
ation domination. This is a combination of the expansion 
deviation acting just like the photons, and a slight prefer- 
ence of the data for additional radiation energy density, 
in accord with previous hints that the number of effec- 
tive neutrino species, N c g, might be greater than the 
standard model value of 3.046. Indeed, the mean value 
of S 2 = 0.026 in this case corresponds to AiV off = 0.31, 
in good agreement with the Planck values of N c s = 3.39. 
In all other parts of the article we keep c s = 1. 
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FIG. 3. Reconstruction of the expansion history deviations 
8(a) from ACDM is shown, with the mean value (solid line) 
and 68% uncertainty band (shaded area). 



IV. PHYSICAL IMPLICATIONS 

This analysis has been model independent, allowing in- 
dividual epochs to float freely without assuming a func- 
tional form. If we do assume a specific model, then con- 
straints will in general be tighter, with each epoch having 
leverage on others through the restricted form. 

Three distinct families of early dark energy might be 
considered: where the early dark energy density rises, 
falls, or stays constant across CMB recombination. These 
were investigated in [3] in terms of the (somewhat moti- 
vated) models of barotropic aether, dark radiation, and 
Doran-Robbers [11] forms, respectively (see [3] for more 
detailed discussion). We compute the constraints on the 
fraction Q e of critical density contributed by early dark 
energy (approximately equivalent to S) within each of 
these models (not using the Si bins), giving the results in 
Table II. (Note that Planck finds fi e < 0.009 at 95% CL 
for the Doran-Robbers model when also including high 
multipole data [12].) 



TABLE I. 95% confidence upper bounds are given for the ex- 
pansion history deviations 5, listed by the bin number and 
midpoint of the log a bins, for cases with different sound 
speeds. 



Figure 3 shows the mean value and 68% uncertainty 
band of the expansion deviations 6(a) given by the Monte 
Carlo reconstruction using the recent CMB data. This 
figure represents the best current model-independent 
knowledge of the early expansion history of our Universe. 
Setting all Si = 0, i.e. ACDM, is consistent with these re- 
sults at the 95% confidence level. The mean value does 
show a very slight preference for a faster expansion rate, 
as in early dark energy or extra relativistic degrees of 
freedom, before recombination. 
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Aether Dark Radiation Doran-Robbers 



0.019 



0.033 



0.012 



TABLE II. The 95% confidence level uncertainties are pre- 
sented for three early dark energy models. For small values, 
Q, e w 8. The Doran-Robbers model has an additional param- 
eter w ; we find w = — 1.49±g;|? (95% CL). 



Two aspects of the models impact their detectability: 
the presence of the expansion history deviation at a sen- 
sitive epoch and its persistence over time, and its clus- 
tering behavior. The common Doran-Robbers form has 
the tightest bounds (despite the extra parameter), due to 
its persistence pre- and post-recombination and its dis- 
tinction from matter clustering since it has c 2 = 1. The 



aether model only begins to deviate around recombina- 
tion, and has c 2 = so there is more covariance with 
the dark matter component. Dark radiation has influ- 
ence only before recombination and its c 2 = 1/3 makes it 
more covariant with the photons (and neutrinos). A key 
conclusion is that early dark energy could in fact be more 
prevalent than apparent from bounds in the literature on 
the Doran-Robbers model. 

Since dark radiation density at early times scales like 
radiation, it acts like the addition of relativistic degrees of 
freedom. Taking into account the definition of extra de- 
grees in terms of the number of effective neutrino species 
iV e fj, the constraint on 5 within the dark radiation model 
translates to [13] 

AJV cff « 7.44<5 darkrad « 7.44fi e /(l - n e ) . (4) 

Thus f2 e < 0.033 for the dark radiation model becomes 
AiVcff < 0.25 at 95% CL. This puts a tighter global 
bound on AiV c ff compared to our model independent 
value from 82 before recombination (AN e g < 0.71 at 95% 
CL). 

Another implication of the expansion history is its re- 
lation to the spacetime itself. The Ricci scalar curvature 
is the central quantity in the Einstein-Hilbert action for 
general relativity, and plays a key role as well in exten- 
sions to gravity such f(R) theories. The curvature his- 
tory of the Universe has been explored from a theoretical 
perspective recently by [14]. Since 
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deviations, whether due to, e.g., some form of dark en- 
ergy or a nonstandard number of relativistic degrees of 
freedom. It also relates directly to the Ricci spacetime 
curvature. 

Adding late time data that helps to constrain Hq or 
f2 m , say, would help break the degeneracy around re- 
combination that led to the loosest, 16% upper bound 
on deviations. However, proper treatment of this would 
require many low redshift bins to reflect the density of 
the data, while our focus here is on the early expansion 
history. 

We regard the model independence of the analysis as 
a signal virtue; however we can also compare the bounds 
for specific early dark energy models. For the barotropic 
aether, dark radiation, and Doran-Robbers models we de- 
rive 95% CL limits of less than 0.019, 0.033, 0.012 in early 
dark energy density 51 e , respectively. We emphasize that 
bounds appear tightest when assuming the conventional 
Doran-Robbers form, and so early dark energy should be 
not be thought ruled out based purely on constraining 
this model. In terms of extra effective neutrino species 
the model independent results imply AN c g < 0.71 at 
95% CL. 

Future CMB data, such as the release of polarization 
data from Planck, ACTpol [15], PolarBear [16], SPTpol 
[17] experiments will enhance our knowledge of the his- 
tory back to z « 10 5 . Exploring the even earlier universe 
will require neutrino, dark matter, or gravitational wave 
astronomy. Late time probes will continue to map the 
last e-fold of cosmic expansion in greater detail. Over just 
a few years cosmological observations have taken us from 
order unity uncertainty (with > 10% in narrow epochs 
around recombination and today) to a few percent level 
knowledge over more than 10 e-folds of cosmic history. 



observational constraints on 8 (and hence w through 
Eq. 3) can be used to cast light on the curvature his- 
tory. 



V. CONCLUSIONS 

We have used the recent advances in CMB data to con- 
strain the fundamental quantity of the expansion history 
of our Universe. The results from the model indepen- 
dent analysis bound deviations from ACDM at 2-16% 
(95% CL), depending on the epoch. This constrains any 
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